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We demonstrate an efficient high-precision optical spectroscopy technique for single trapped ions 
with non-closed transitions. In a double-shelving technique, the absorption of a single photon is 
first amplified to several phonons of a normal motional mode shared with a co-trapped cooling ion 
of a different species, before being further amplified to thousands of fluorescence photons emitted 
by the cooling ion using the standard electron shelving technique. We employ this extension of 
the photon recoil spectroscopy technique to perform the first high precision absolute frequency 
measurement of the ^ 03/2 ^ ^Pi /2 transition in ^°Ca^, resulting in a transition frequency of 
f — 346 000 234 867(96) kHz. Furthermore, we determine the isotope shift of this transition and the 
^Si /2 ^ ^Pi /2 transition for ^^Ca^, ^^Ca^ and ^®Ca^ ions relative to ^°Ca^ with an accuracy 
below 100 kHz. Improved field and mass shift constants of these transitions as well as changes in 
mean square nuclear charge radii are extracted from this high resolution data. 

PACS numbers: 42.50.-p, 03.75.Be, 32.10.Bi, 21.10.Ft 


Single or few laser-cooled ions trapped in Paul traps al¬ 
low spectroscopy with high resolution and accuracy due 
to the near perfect realization of an unperturbed tran¬ 
sition frequency as demonstrated in ion-based optical 
clocks [1]. For long-lived excited states the electron shelv¬ 
ing technique provides near unity detection efficiency [ 2 ] . 
Advanced laser fluorescence techniques are well-suited for 
closed or nearly-closed broad, dip ole-allowed transitions 
(am- However, they provide insufficient signal to effi¬ 
ciently study non-closed transitions where typically only 
one photon can be absorbed before the ion decays to a 
state which is not addressed by the spectroscopy laser. 
As a result, the fluorescence rate is low even when pump- 
repump techniques are employed, resulting in long av¬ 
eraging times 0 E]. Recently, new spectroscopy tech¬ 
niques with a high sensitivity to the absorption of single 
or few photons have been developed employing the signal 
from photon recoil mm or decoherence of an electronic 
superposition state [9]- An evaluation of the dominant 
effects shifting the resonance from its unperturbed fre¬ 
quency provides accuracies below 100 kHz m [ 8 ]. Up to 
now, these absolute frequency measurements have been 
demonstrated on nearly-closed transitions. 

Here, we demonstrate an extension of the photon recoil 
spectroscopy (PRS) technique [ 8 ] to perform an abso¬ 
lute frequency measurement of the non-closed ^ 03/2 ^ 
^Pi /2 transition in ^^Ca+ with an accuracy of below 
100 kHz. Using this extension together with the original 
PRS scheme, we performed isotope shift measurements 
of the ^D 3/2 ^ ^Pi /2 aiid ^Si /2 ^ ^Pi /2 transitions 
for the most abundant even Ca+ isotopes with an uncer¬ 
tainty of below 100 kHz. These measurements improve 
the accuracy of the isotope shifts of the Ca+ ions by up 
to two orders of magnitude compared to previous mea¬ 
surements nnmi]. 

Prom a multi-dimensional King plot analysis na [13], 


seeded by the changes in mean square nuclear charge 
radii obtained from muonic atom spectroscopy and elec¬ 
tron scattering, we derive improved values for the field 
and mass shift constants as well as the change in the 
mean square nuclear charge radius. The high accuracy 
of these measurements serves as a benchmark for ab initio 
atomic structure calculations and enables improved cal¬ 
ibration for collinear laser spectroscopy mm, where 
typically uncertainties of a few MHz are achieved. These 
measurements are relevant for the investigation of nu¬ 
clear properties of unstable calcium isotopes usms]. 
Calcium is of particular interest since its isotopes include 
two doubly magic nuclei at A = 40 and 48. In the shell 
model, the additional 8 neutrons form a complete I/ 7/2 
shell, however there is little change in the mean-square 
proton radius (charge radius). Understanding these pe¬ 
culiar changes in the calcium charge radius is a challenge 
for nuclear theory. 

We use the experimental setup for PRS with ^^Mg+ as 
the cooling ion species described in detail in [HiiiaiiH]. 
In brief, we prepare two-ion crystals in a linear Paul trap 
consisting of ^^Mg+ and a selected Ca+ isotope as the 
spectroscopy ion. Isotope-selective two-photon ionization 
is used to load the different Ca+ isotopes: the ^Si ^ 
^Pi transition of the neutral species of the target isotope 
is selectively excited using a narrow-linewidth laser m, 
followed by ionization with a second laser. After load¬ 
ing, the isotopes are identified by the motional normal 
mode frequencies uji^o (h iii phase, o: out of phase) of 
the two-ion crystal in the axial trapping potential, given 
by [inilll], 

(1 + M T \/l -M + M^) , (1) 

where ujz is the axial normal mode frequency of a sin¬ 
gle ^^Mg+, /i = '^Mg/'^Ca? and rriMg and mca are the 
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FIG. 1. Experimental sequence for the spectroscopy on (a) 
the ^Si/2 —^ ^Pi/2 transition and. (b) the ^03^2 —^ "Pl/2 
transition. A detailed description is given in the text, (c) 
Simplified level scheme of Ca^ ions 


masses of the two ions in the trap. Following the tech¬ 
nique described in [ 20 ], the normal mode frequencies were 
measured and fits of the in-phase normal mode frequen¬ 
cies lead to 1.922(1) MHz, 1.889(1) MHz, 1.857(1) MHz 
and 1.795(1) MHz for a two-ion crystal containing ^^Ca+, 
^^Ca+, "^^Ca^ and ^^Ca+, respectively. The measured 
frequencies are in agreement with the frequencies calcu¬ 
lated from 0 for the measured axial normal mode fre¬ 
quency of a single ^^Mg+ of 2.221(1) MHz. 

Isotope shift measurements for the ^Si /2 ^ ^Pi /2 tran¬ 
sition of the different Ca+ isotopes were performed using 
the original PRS scheme [8]. The experimental sequence 
of the technique is displayed in Fig. 0 )- After Doppler 
cooling, we further cool the axial modes of the two-ion 
crystal to the motional ground state via sideband cool¬ 
ing on ^^Mg+ [18]. Starting from the motional ground 
state of one of the axial normal modes of the two-ion 
crystal, the spectroscopy laser with a wavelength of 397 
nm is applied in 125 ns short pulses. The laser pulse 
repetition is matched to the motion of the two-ion crys¬ 
tal, thus enhancing the sensitivity through resonant driv¬ 
ing and enabling the detection of around 10 absorbed 
photons with a signal-to-noise ratio of 1 [8]. A repump 
laser at 866 nm resonant with the ^ 03/2 ^ ^Pi /2 tran¬ 
sition is applied interleaved with the spectroscopy laser 
pulses to prevent populating the metastable ^D 3 / 2 -state. 
These spectroscopy/repump cycles are repeated 70 times, 
obtaining the maximum motional excitation detectable 
without saturating the signal. The residual ground state 
population is measured using a red sideband stimulated 
Raman adiabatic passage (STIRAP) pulse sequence cou¬ 
pling the two ground state hyperfine levels in ^^Mg+ 
|22| . This maps the motional excitation efficiently into 
electronic excitation, detectable via discrimination of the 
25Mg+ hyperfine states using electron shelving with the 
TT-detection technique [53]. 

We extend the original PRS scheme to measure the ab¬ 
solute frequency and the isotopic shifts of the open ^D 3/2 
^ ^Pi /2 transition in Ca+ as shown in Fig. mb). The 
sequence starts by initializing the Ca+ ion in me ^ 03/2 
state via optical pumping on the ^Si /2 ^ ^Pi /2 tran¬ 


sition and subsequent sideband cooling on ^^Mg+. Ab¬ 
sorption of a photon on the ^D 3/2 ^ ^Pi /2 spectroscopy 
transition in Ca+ results in a decay to the ^Si /2 ground 
state with 94 % probability [24], where the ion is lost from 
the normal spectroscopy cycle. Starting from the ground 
state of motion of one of the ion’s normal modes, absorp¬ 
tion of photons from a pulsed laser tuned to resonance 
with the ^Si /2 ^ ^Pi /2 transition in Ca+ efficiently ex¬ 
cites this motional mode analogously to the original PRS 
technique [ 8 ]. This corresponds to a first shelving step, 
since the excitation of phonons is conditional upon the 
absorption of a photon on the spectroscopy transition. 
Due to the non-vanishing probability of the ion to de¬ 
cay to the ^D 3/2 state the ion is effectively reinitialized 
in this state after 70 excitation pulses. Repeating the 
spectroscopy/motional excitation cycle three times fur¬ 
ther enhance the motion of the two-ion crystal to obtain 
the maximum motional excitation without reaching sat¬ 
uration. The motional excitation is mapped onto the 
25Mg+ ion and detected as before. 

Commercially available external cavity diode lasers are 
used to excite the transitions of the different Ca+ iso¬ 
topes, where the laser coupling the ^Si /2 and the ^Pi /2 
state is frequency doubled to 397 nm. Both laser beams 
used in the two spectroscopy techniques are aligned span¬ 
ning a 45° angle with an external magnetic field required 
for efficient cooling and detection. The polarization is 
aligned perpendicular to the magnetic field to couple all 
magnetic substates of the involved transitions. The fre¬ 
quencies of both lasers are stabilized in the infrared to 
the frequency doubled output of a fiber-based optical fre¬ 
quency comb using an electronic feedback loop. A rela¬ 
tive frequency stability of better than 10 “^^ is achieved 
with respect to the comb for averaging times of a second. 
The comb is stabilized to the 10 MHz reference signal ob¬ 
tained from a hydrogen maser, which is referenced to the 
SI second at the German Metrology Institute (PTB). 
The lasers are steered to the center frequencies of the 
spectroscopy transitions using acousto-optical modula¬ 
tors in double-pass configuration. We apply the two- 
point sampling technique [ 8 ], where the population dif¬ 
ference at half maximum below and above the resonance 
line is measured and the difference is used to correct the 
center of the two probing frequencies for the next scan. 
The population difference in combination with the ex¬ 
perimentally determined slope at these positions is used 
in post evaluation to derive the center frequency of the 
atomic transition. 

Following [ 8 ] we evaluated the dominant effects which 
could result in a shift of the observed transition frequency 
from its unperturbed value. External magnetic fields in 
combination with polarization asymmetries can cause an 
asymmetric line shape shifting the center frequency of 
the transition. Furthermore, the electrical field from the 
spectroscopy laser light may shift the frequency by means 
of the ac stark shift, even though small intensities on the 
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FIG. 2. Overlapped Allan deviation cry{r) of the frequency 
measurement of the ^D 3/2 ^ ^Pi /2 transition of ^°Ca^. The 
inset shows a resonance scan of the transition, where the red 
line corresponds to density matrix simulations adjusted in am¬ 
plitude and offset to account for experimental imperfections. 


order of a few % of the saturation intensity were em¬ 
ployed for probing. Therefore, we measured the tran¬ 
sition frequency for different magnetic fields and differ¬ 
ent spectroscopy laser intensities and extracted the cor¬ 
responding shifts A/zeeman = 54 ± 44 kHz and A/stark = 
—70 ± 83 kHz for the experimentally used parameters. 
The given uncertainties correspond to the 68.3 % predic¬ 
tion bound of the fit to the measured data and are limited 
by statistics. Due to the first shelving step in the newly 
developed technique, the dominant shift in the PRS tech¬ 
nique, the so-called lineshape shift arising from probe- 
frequency-dependent Doppler cooling and heating effects 
during motional excitation, is expected to be small. This 
was confirmed using density matrix simulations, where 
the shift was determined to be well below 1 kHz. The 
absolute frequency of the ^D 3/2 ^ ^Pi /2 transition was 
derived as the weighted average of 6 frequency measure¬ 
ments and is given by / = 346 000 234 867(96) kHz, where 
94 kHz of the uncertainty results from systematic and 21 
kHz from statistical uncertainties. 

We used the original and extended PRS technique to per¬ 
form isotope shift measurements for the ^Si /2 ^ ^Pi /2 
and ^D 3/2 ^ ^Pi /2 transitions, respectively. The mea¬ 
surements of each isotope and the reference ^^Ca+ were 
performed interleaved. We first measured the "^^Ca^ 
transition frequency. Then, we loaded a new two-ion 
crystal containing the investigated isotope of Ca+ and 
measured its transition frequency. After this, we re¬ 
peated the measurement of the reference ion ^^Ca+ to 
ensure the absence of frequency drifts. 

We investigated the statistical properties of all frequency 
measurements by plotting the overlapped Allan devia¬ 


tion cry{r). No indication of drifts or other frequency 
fluctuations were observed for the performed measure¬ 
ments. As an example, the overlapped Allan deviation 
for the measurement of the ^D 3/2 ^ ^Pi /2 transition 
of ^0Ca+ is shown in Fig. After averaging, the mea¬ 
sured frequencies were used to derive the isotope shift 
for each measurement. For the ^Si /2 ^ ^Pi /2 transi¬ 
tion the measured isotope shifts were corrected by the 
change of the probe frequency-dependent lineshape shift, 
resulting from the different mass of the isotope and the 
corresponding change of the normal mode frequency of 
the two-ion crystal. Using the simple model derived in 
[8], the calculated correction is smaller than 2 kHz per 
mass unit corresponding to a maximum correction of 10 
kHz. All other systematic effects are common for the 
measurements performed on the same day day and are 
therefore removed. The final values with the cor¬ 

responding accuracies are displayed in Table [T[ 

Isotope shifts in atomic transition frequencies originate 
from two effects: a change in the size of the nucleus and 
thus its interaction with the electrons (field shift), and a 
change in the recoil of the nucleus (mass shift). The mass 
shift is more important for light elements, while for heavy 
elements the field shift dominates [25]. The difference in 
the transition frequency, , between isotopes with 

mass rriA and rriA' can be expressed as 

where kus is the mass shift constant, F is the field shift 
constant, and is the change in the mean-square 

nuclear charge radius. The latter is a common parameter 
when comparing two different transitions and can be 
eliminated in a King plot analysis as shown in Fig. 
for the two transitions considered here. Each axis 
shows the modified isotope shift , 

where = (1/m^ — l/mA ')~^for one of the two 

transitions. A straight line fit to the three data points 
provides linear combinations of the field and mass shift 
constants for the two transitions. An important result 
from this fit is that there is no evidence for a deviation 
from a straight line, confirming that is a good 
parametrization of the isotope shift even at the high 
experimental accuracy of the measurements presented 
here. 

A comparison of the high resolution results with previous 
experimental data based on collinear laser spectroscopy 


TABLE I. Measured isotope shifts of the two transitions ref¬ 
erenced to ^°Ca^. 


A 

54’^ (MHz) 

(MHz) 

42 

425.706 (94) 

-2349.974 (99) 

44 

849.534 (74) 

-4498.883 (80) 

48 

1705.389 (60) 

-8297.769 (81) 
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FIG. 3. Two-dimensional King plot showing modified isotope 
shift of the 866 and 397 lines. Red squares: previous ex¬ 
perimental data from US] and m-, blue circles: this work. 
The insets show the relevant ranges enlarged by a factor of 
approximately 30 to illustrate the quality of the fit. 


nnuni shows systematic deviations, which can be used 
to calibrate experimental parameters of this technique. 
Following Ref. na we performed a three-dimensional 
King plot analysis to extract the fitting parameters 
kus and F for the two transitions. Two dimensions 
are those shown in Fig. [3 In the third dimension we 
plot the modified change in mean-square nuclear charge 

radius using the previous values of 5(r‘^) 

from [26], which are based on muonic atom spectroscopy 
and electron scattering. The three-dimensional King 
plot constrains the mass and field-shift constants, and 
under the assumption that is correct (i.e. the three 
data points are connected by a straight line) can also 
be used to extract improved values of To find 

the parameter estimates and their uncertainties an 
acceptance/rejection Monte Carlo method was used to 
generate samples consistent with the measured values 
and associated uncertainties. The measurement distri¬ 
butions were assumed to be independent uncorrelated 
normals. The likelihoods of three randomly generated 
points, constrained to be collinear, were used as the 
acceptance criterion in the algorithm. The extracted 
parameters are shown in Table [TT] 

The extracted field-shift and mass-shift constants 
pose a strong challenge for many-body atomic theory 
(fourth column of Table [^, where the mass-shift in 
particular has proven very difficult to calculate even in 
the “easy” case of single-valence-electron ions [izlEHI. A 
comparison to the experimental field and mass shift con¬ 
stants given in nnmi] proves difficult since the derived 
uncertainties depend strongly on the analysis technique 
and input parameters for Evaluating the field and 


TABLE II. Parameters of three-dimensional King plot seeded 
with values of taken from |26| . The units for the 

field Fi and mass ki shift constants and the changes in mean 
square nuclear charge radii are MHz-fm“^, GHz-amu 

and fm^, respectively. For comparison the second column for 
the previous data shows results for the analysis using isotope 
shift data taken from m and in analyzed with the methods 
used in this work. 


Par am. 

Previous 

This work 

Theory 

F 397 

-283(6Q 

-281(34) 

-281.8(7.0) 

-285(^ 




-28iT 

k397 

4 ^ 

0 

or 

bo 

406.4(2.8) 

408.73(40) 

35^ 

427"' 

Fsgg 

79(4:(] 

80(13) 

87.7(2.2) 

8^ 

92 

-220 (F 

ksGG 

-1989.8(4;|^ 

-1990.9(1.4) 

-1990.05(13) 

-218f'^ 

./ 2\42,40 

(5(r ) 

0.210(7) 

0.210(7) 

0.2160(49) 


r/ 2 \ 44,40 

(5(r ) 

0.290(9) 

0.290(9) 

0.2824(65) 


c/ 2\48,40 

S{r ) 

-0.005(6) 

-0.005(6) 

-0.0045(60) 



^ Martensson-Pendrill et al.^ 1992 |10| . 

^ Safronova and Johnson, 2001 |27| . 

Nortershauser et a/., 1998 [77] . 

^ This work, based on the methods in EH). 


mass shift constant from isotope shifts given in HQiin] 
using the analysis and input parameters employed here, 
indicates an order of magnitude improvement in these 
values based on the high resolution data presented here. 
Furthermore, the analysis leads to improved values of 
changes in mean-square nuclear charge radii which are 
in good agreement with the values derived from isotope 
shift measurements in neutral Ca m- 
In conclusion, we presented the first high precision 
absolute frequency measurement of the ^ 03/2 —> ^Pi /2 
transition in Ca+ using an extension of photon recoil 
spectroscopy. We measured isotope shifts of this and 
the ^Si /2 ^ ^Pi /2 transition of the most abundant 
stable, even isotopes in Ca+ using this highly sensitive 
technique and improved the accuracy of the shifts by up 
to two orders of magnitude. From these measurements, 
improved values of the field and mass shift constants 
as well as changes in mean square nuclear charge radii 
are extracted. These results may serve to validate cal¬ 
culations and determine optimal models used in isotope 
shift calculations. Furthermore, the precise knowledge 
of the isotopic shifts can be used to improve the 
calibration of collinear laser spectroscopy experiments 
na ng and improve the analysis of quasar absorption 
spectra aimed at searches for variations of fundamental 
constants [29|. This work demonstrates that photon 
recoil spectroscopy is a versatile technique, suitable for 
a wide range of high precision spectroscopic applications. 

We acknowledge the support of DFG through QUEST 






























5 


and grant SCHM2678/3-1. This work was financially 
supported by the State of Lower-Saxony, Hannover, Ger¬ 
many. Y. W. acknowledges support from IGSM. We 
thank R. Blatt for generous loan of equipment, W. 
Nortershauser and C. Elster for stimulating discussions 
concerning the data evaluation, and J. B. Wiibbena, 1. D. 
Leroux and N. Scharnhorst for technical assistance with 
the absolute frequency measurement. 


[1] A. D. Ludlow, M. M. Boyd, J. Ye, E. Peik, and P. O. 
Schmidt, arXiv:1407.3493 (2014) 

[2] H. G. Dehmelt, Bull. Am. Phys. Soc. 20, 60 (1975). 

[3] A. Gardner, K. Sheridan, W. Groom, N. Seymour-Smith, 
and M. Keller, Appl. Phys. B 117, 755 (2014) 

[4] M. Herrmann, V. Batteiger, S. Kniinz, G. Saathoff, 

T. Udem, and T. Hansch, Phys. Rev. Lett. 102, 013006 
(2009) 

[5] T. Pruttivarasin, M. Ramm, and H. Haffner, J. Phys. B 
47, 135002 (2014) 

[6] J. B. Wiibbena, Controlling motion in quantum logic clocks , 
Ph.D. thesis, Gottfried Wilhelm Leibniz Universitat 
Hannover (2014). 

[7] C. Hempel, B. P. Lanyon, P. Jurcevic, R. Gerritsma, 

R. Blatt, and C. F. Roos, Nat. Photonics 7, 630-633 
(2013) 

[8] Y. Wan, F. Gebert, J. B. Wiibbena, N. Scharnhorst, 

S. Amairi, 1. D. Leroux, B. Hemmerling, N. Lorch, 

K. Hammerer, and P. O. Schmidt, Nat. Commun. 5, 
4096 (2014) 

[9] G. Clos, M. Enderlein, U. Warring, T. Schaetz, and 
D. Leibfried, Phys. Rev. Lett. 112, 113003 (2014) 

[10] A.-M. Martensson-Pendrill, A. Ynnerman, H. Warston, 

L. Vermeeren, R. E. Silverans, A. Klein, R. Neugart, 

C. Schulz, P. Lievens, and 1. Collaboration, Phys. Rev. 

A 45, 4675 (1992) 


[11] W. Nortershauser, K. Blaum, K. Icker, P. Muller, 
A. Schmitt, K. Wendt, and B. Wiche, Eur. Phys. J. 
D 2, 33 (1998) 

[12] C. W. P. Palmer, P. E. G. Baird, S. A. Blundell, J. R. 
Brandenberger, C. J. Foot, D. N. Stacey, and G. K. 
Woodgate, J. Phys. B 17, 2197 (1984) 

[13] 1. Angeli and K. Marinova, Atomic Data and Nuclear 
Data Tables 99, 69 (2013) 

[14] R. E. Silverans, G. Borghs, P. De Bisschop, and 

M. Van Hove, Hyperfine Interact. 24, 181-201 (1985) 

[15] R. F. Garcia Ruiz, “Ground-State Electromagnetic Mo¬ 
ments of Calcium Isotopes,” (2015), submitted. 

[16] C. Gorges, (2015), priv. comm. 

[17] B. Hemmerling, F. Gebert, Y. Wan, D. Nigg, 1. V. Sher- 
stov, and P. O. Schmidt, Appl. Phys. B 104, 583 (2011) 

[18] Y. Wan, F. Gebert, F. Wolf, and P. O. Schmidt, preprint 

arXiv:1501.06987 (2015), accepted for publication in 

Phys. Rev. A. 

[19] D. Lucas, A. Ramos, J. Home, M. McDonnell, 
S. Nakayama, J.-P. Stacey, S. Webster, D. Stacey, and 
A. Steane, Phys. Rev. A 69, 012711 (2004) 

[20] M. Drewsen, A. Mortensen, R. Martinussen, P. Staanum, 
and J. S0rensen, Phys. Rev. Lett. 93, 243201 (2004) 

[21] J. B. Wiibbena, S. Amairi, O. Mandel, and P. O. 
Schmidt, Phys. Rev. A 85, 043412 (2012) 

[22] F. Gebert, (2015), in preparation. 

[23] B. Hemmerling, F. Gebert, Y. Wan, and P. O. Schmidt, 
New J. Phys. 14, 023043 (2012) 

[24] M. Ramm, T. Pruttivarasin, M. Kokish, 1. Talukdar, and 
H. Haffner, Phys. Rev. Lett. Ill, 023004 (2013) 

[25] G. Breit, Rev. Mod. Phys. 30, 507 (1958) 

[26] H. D. Wohlfahrt, E. B. Shera, M. V. Hoehn, Y. Yamazaki, 
G. Fricke, and R. M. Steffen, Phys. Rev. C 73B, 131 
(1978). 

[27] M. S. Safronova and W. R. Johnson, Phys. Rev. A 64, 
052501 (2001). 

[28] J. C. Berengut, V. A. Dzuba, and V. V. Flambaum, 
Phys. Rev. A 68, 022502 (2003). 

[29] M. T. Murphy and J. C. Berengut, MNRAS 438, 388 
(2014) 




